SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWlwl  Dim  Entered; 


REPORT  DOCUMENTATION  PAGE 


REPORT  NUMBER 

NRL  Memorandum  Report  4389 


«.  TITLE  (and  Submit) 

NUMERICAL  SIMULATIONS  OF  EQUATORIAL 
SPREAD  F  USING  ALTAIR  INCOHERENT  BACKSCATTER 
MEASURED  ELECTRON  DENSITY  PROFILES  FROM  THE 
17  JULY  1979  DNA  KWAJALEIN  CAMPAIGN 


7.  AUTMORCJ 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


7  RECIPIENT'S  CATALOG  NUMBER 


S.  TYPE  OF  REPORT  «  PERIOD  COVERED 

Interim  report  on  a  continuing 
NRL  problem. 


S.  PERFORMING  ORG.  REPORT  NUMBER 


t.  CONTRACT  OR  GRANT  HUMBERT*; 


S.  T.  Zalesak  and  S.  L.  Ossakow 


S.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Research  Laboratory 
Washington,  DC  20375 


<0.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  S  WORK  UNIT  NUMBERS 


67-0889-0-0 


<1.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


Defense  Nuclear  Agency 
Washington,  DC  20305 


IZ.  REPORT  DATE 

December  3, 1980 


I S.  NUMBER  OF  PAGES 
28 


MONITORING  AGENCY  NAME  A  AODRESSf/f  different  from  Confroflln*  Office J  IS.  SECURITY  CLASS,  fof  tbit  report; 

UNCLASSIFIED 


•«.  DISTRIBUTION  STATEMENT  fof  tbit  Report; 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  tho  mbetrect  motored  In  Block  20,  If  different  from  Report) 


IB.  SUPPLEMENTARY  NOTES 

This  research  was  sponsored  by  the  Defense  Nuclear  Agency  under  subtask  S99QAXHC041, 
work  unit  21,  and  work  unit  title  “Plasma  Structure  Evolution.” 


19.  KEY  WORDS  fContJnu*  on  reverie  mldo  if  nocooeory  end  Identify  by  block  number) 

Equatorial  Spread  F 
Numerical  simulations 
Collisional  Rayleigh-Taylor  mechanism 
ALTAIR  electron  density  profiles 
Spread  F  bubbles 


20.  ABSTRACT  fContlnue  on  reeeree  tide  II  neceeeery  and  Identify  by  block  numbtr) 

Using  electron  density  profiles  measured  by  the  ALTAIR  radar  on  17  July  1979,  we  have  per¬ 
formed  global  large  scale  size  numerical  simulations  of  the  nonlinear  evolution  of  the  collisional 
Rayleigh-Taylor  instability  in  the  equatorial  ionosphere.  The  ALTAIR  profiles  were  taken  approxi¬ 
mately  four  and  one-half  hours  prior  to  the  launch  of  a  rocket  equipped  with  plasma  diagnostics 
instrumentation,  and  on  the  order  of  one  hour  prior  to  the  onset  of  equatorial  spread  F.  Using  5% 
amplitude  sinusoidal  initial  perturbations  in  our  numerical  simulations,  we  find  fully  developed 
spread  F  bubbles  (plumes)  on  time  scales  of  approximately  one-half  hour  for  both  small  (8  km)  and 
luge  (200  km)  horizontal  scale  lengths.  ■*- 


EDITION  OF  I  NOV  ••  IS  OBSOLETE 
i/H  01 02-LF-01 4-6601 


_-rce 


L  J  k  t> 


security  classification  OF  THIS  RAOl 


~~  v> 


CONTENTS 

INTRODUCTION  . 

THEORY . 

NUMERICAL  SIMULATION  RESULTS  AND  DISCUSSION 

SUMMARY . 

ACKNOWLEDGMENTS . 

REFERENCES  . 


Accession  For  ] 

MTIS  GRA&I 

MDC  TAB 

Unannounced 

Justification 

\ 

i 

NUMERICAL  SIMULATIONS  OF  EQUATORIAL  SPREAD  F  USING  ALTAIR 
INCOHERENT  BACKSCATTER  MEASURED  ELECTRON  DENSITY  PROFILES 
FROM  THE  17  JULY  1979  DNA  KWAJALEIN  CAMPAIGN 

I.  INTRODUCTION 

During  the  17  July  1979  DNA  PLUMEX  I  experiment,  coordinated  measurements 
of  equatorial  spread  F  (ESF)  by  radar  backscatter  using  the  ALTAIR  radar,  by 
plasma  diagnostics  using  in  situ  rocket  probes,  and  by  scintillation  measure¬ 
ments  using  the  DNA  Wideband  satellite,  were  carried  out.  Prior  to  the  rocket 
launches,  however,  ALTAIR  was  used  to  monitor  the  equatorial  environment  into 
which  the  rockets  would  be  launched.  In  particular,  off-perpendicular 
"incoherent"  scans  were  made,  giving  direct  profiles  of  electron  density 
versus  altitude,  as  well  as  coherent  backscatter  measurements.  Using  a 
background  laminar  electron  density  profile  of  this  type,  taken  before  the 
onset  of  spread  F,  as  the  ambient  laminar  ionospheric  profile,  we  have 
performed  numerical  simulations  of  the  nonlinear  evolution  of  the  collisonal 
Rayleigh- Taylor  instability.  This  phenomenon  is  believed  to  be  the  cause1-5 
on  the  large  scale  irregularities  during  ESF. 

Our  results  indicate  fully  developed  spread  F  plumes  (bubbles)  on  time 
scales  of  approximately  half  an  hour,  using  initial  sinusoidal  perturbations 
of  5%  amplitude  and  two  different  horizontal  wavelengths  (8  km  and  200  km). 
ALTAIR  radar  and  ionosonde  data  indicate  fully  developed  spread  F  approxi¬ 
mately  one  hour  after  the  aforementioned  profile  was  taken6.  The  difference 
in  times  could  easily  be  accounted  for  simply  by  assuming  a  smaller  initial 
perturbation,  i.e.,  the  present  measurements  do  not  provide  us  with  all  the 
pertinent  initial  conditions  prior  to  ESF  onset.  Other  factors,  such  as 
significant  differences  between  local  and  magnetic  field  line  integrated 
Pedersen  conductivities,  the  shorting  effects  of  background  E  region  con¬ 
ductivities,  and  the  effect  of  inertial  terms  in  the  ion  momentum  equation 

(ignored  in  our  simulations)  could  also  delay  the  progress  of  the  Instability. 
Manuscript  submitted  October  6, 1980. 


Notwithstanding  these  effects,  our  aim  was  to  show  that  these  types  of  global 
large  scale  numerical  simulations  of  the  collisional  Rayleigh-Taylor  mechanism 
are  consistent  with  the  onset  time  and  morphology  of  ESF  during  the  17  July 
1979  campaign.  Also,  we  find  that  the  larger  horizontal  scale  length  bubbles 
consist  of  plasma  which  originated  at  much  lower  altitudes  than  that  of 
their  smaller  horizontal  scale  length  counterparts,  resulting  in  much  more 
depleted  bubbles  in  the  larger  horizontal  scale  length  case-  This  last  effect 
is  explained  using  scaling  arguments  similar  to  those  which  apply  to  the 
fringe  field  at  the  edge  of  a  parallel  plate  capacitor.5 

In  Section  11  we  briefly  review  the  relevant  theory  and  equations 
used  in  the  numerical  simulations;  in  Section  III  we  present  and  analyze 
the  numerical  results;  and  in  Section  IV  a  summary  is  presented. 

II.  THEORY 

For  a  complete  description  of  our  theoretical  and  numerical  model  of 
the  equatorial  spread  F  phenomenon  see  reference  4.  Briefly,  we  assume 
that  the  magnetic  field  line  integrated  Pedersen  conductivity  in  the 
equatorial  ionosphere  is  dominated  by,  and  is  therefore  proportional  to, 
the  local  (equatorial)  Pedersen  conductivity,  which  in  turn  is  proportional 
to  the  product  of  the  local  ion-neutral  collision  frequency  and  the  local 
electron  density.  This  enables  us  to  carry  out  our  simulations  in  a  two 
dimensional  (x,y)  coordinate  system.  The  constant  magnetic  field  B  is 
aligned  along  the  z  axis  (pointing  north) .  Gravity  is  directed  along  the 
negative  y  axis.  n(y),  VR(y),  and  V^n(y)  are  the  ambient  electron  density, 
recombination  coefficient,  and  ion-neutral  collision  frequency  respectively. 
Magnetic  field  lines  are  assumed  to  terminate  at  both  ends  in  an  electrically 


insulated  medium  (currents  must  close  in  the  two  dimensional  plane,  not 
in  some  distant  E  region). 

Following  reference  4,  we  describe  the  system  with  the  two-fluid  plasma 
continuity  and  momentum  equations: 


dn 

+  V  •  (n  v  )  =  -  V  n 
at  a-a  R  a 


(1) 


(2) 


where  the  subscript  a  denotes  the  species  (i  for  ions,  e  for  electrons), 

n  is  the  species  number  density,  v  is  velocity,  VR  is  the  recombination 

coefficient,  E  is  the  electric  field,  £  is  the  gravitational  acceleration, 

q  is  the  species  charge,  V  is  the  species  collision  frequency  with  the 

cm 

neutral  atmosphere,  is  the  neutral  wind  velocity,  c  is  the  speed  of 
light,  and  m  is  the  species  mass. 

Note  that,  in  contrast  to  reference  4,  we  have  dropped  the  term  + 

v  n 

R  ao  from  (1).  This  is  the  equivalent  of  dropping  the  assumption  of  the 
existence  of  an  ionization  source  given  by  that  term.  This  ionization 
source  was  such  that  the  ambient  ionization  profiles  n0„(y)  was  an  equili¬ 
brium  profile  (dna0/5t  *  0).  Our  present  model  therefore  has  instead 
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Hence,  when  normalized  resulcs  n^(x,y)/nao(y)  are  later  presented,  it  should 
be  understood  that  both  the  numerator  and  denominator  are  time  dependent. 

Figure  1  shows  the  recombination  rate  and  ion-neutral  collision 

frequency  used  in  our  simulations.  It  shall  be  seen  presently  that  Ven 

need  not  be  specified  as  long  as  it  is  much  smaller  chan  the  electron  gyro- 

frequency  For  details  on  the  form  of  V^n  and  as  depicted  in  Figure 

1,  see  reference  4.  If  we  neglect  the  Inertial  terms  (the  left-hand  side) 

of  (2)  by  assuming  the  Inertial  time  scales  are  much  larger  than  either 

the  gyro-periods  or  the  mean  time  between  collisions,  then  the  equation  can 

be  inverted  to  give  an  algebraic  expressions  for  v^.  In  two-dimensional 

(x,y)  geometry  with  B  along  the  z  axis,  the  solution  is  for  our  problem, 

with  m  «  m,,  VJ  /fl_,  <<  1,  V  /fl  =0  (where  «  eB/m  c) ,  and  U  *  0, 
e  i  in  i  ene  a  a  -n 

B 

^e  "f  £  X  *•  *  “TIT  (4) 


Vi-  (§-  +  f  E)  xi+^L(|-  +  f  E)  (5) 

i  x  i 

We  now  make  the  electrostatic  approximate, 

E  ■  V  (6) 

A  9 

A  A 

where  £  x(5/5x)  +  y(5/Sy),  and  the  quasi-neutrality  approximation 

n  »  n,  £  n.  We  then  have 
e  i 


V  •  J  -  0  (7) 

i  a  en  (v,  -  v  )  (8) 
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ION-NEUTRAL  COLLISION  FREQUENCY,  uin  (sec 


Substituting  (4)  and  (5)  into  (8)  and  evaluating  (7),  we  have  for  the 
electrostatic  potential: 


’j.  •  <Vln”  'x  ♦> 


i  6  /  v  B  dn 

—  g  (V.  n) - g  "r- 

e  BSy  in  c  °  dx 


(9) 


As  in  reference  4  we  set  $  =  <t>0  +  <t»  1  where  V^o  -  -  (tn^g/e)  y-  Since  V^2  <t>Q  *  0, 


our  final  potential  equation  becomes 


V  •  (v ,  n  V  4>, )  *  -  ^  I51 
j.  in  a  1  c  dx 


(10) 


The  effect  of  <j>o  is  merely  to  superimpose  a  bulk  westward  plasma  velocity 
g/fi^  on  the  electron  velocity  field  determined  from  (Jij,  without  affecting 
the  morphology  of  the  developing  structures.  Hence,  we  ignore  this 
motion.  In  addition,  we  have  ignored  any  upward  or  downward  bulk  motion 
of  the  ionospheric  plasma,  i.e.,  any  ambient  eastward  or  westward  electric 
field.  This  is  done  because  we  have  chosen  a  starting  time  for  our 
simulation,  which  while  prior  to  ESF  onset,  which  coincides6  with 
essentially  no  upward  or  downward  bulk  motion  of  the  ionospheric  plasma. 

Our  assumption  of  quasi-neutrality  has  made  one  of  our  two  continuity 
equations  (1)  redundant.  We  therefore  choose  the  electron  equation  for  its 
simplicity : 


III.  NUMERICAL  SIMULATION  RESULTS  AND  DISCUSSION 
Equations  (10)  and  (11),  together  with  appropriate  boundary  conditions, 
constitute  the  nonlinear  system  of  equations  we  shall  solve  numerically. 

The  numerical  calculations  to  be  presented  were  performed  on  a  two-dimensional 
cartesian  (x,y)  mesh  using  42  points  in  the  x  (east-west)  direction,  and 
142  points  in  the  y  (vertical)  direction.  The  (uniform)  grid  spacing  was 
2  km  in  the  y  direction  for  all  calculations.  The  grid  spacing  in  the  x 
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direction  was  200m  in  the  "small"  horizontal  scale  length  cases  and  5  km  in 
the  "large"  cases.  The  bottom  of  the  grid  corresponds  to  282  km  altitude 
and  the  top  of  the  grid  to  564  km  altitude  in  all  simulations.  Periodic 
boundary  conditions  were  imposed  on  both  n  and  (Jij  in  the  x-direction.  In  the 
y  direction  transmittive  boundary  conditions  were  imposed  on  n(£jn/6y  =  0)  and 
Neumann  (d4>1/dy  =  0)  boundary  conditions  were  imposed  on  4>^. 

Three  kinds  of  plots  will  be  presented:  (1)  contours  of  constant 
n(x,y,t);  (2)  contours  of  constant  n(x,y,t)/no(y,t) ;  and  (3)  contours  of  con¬ 
stant  electrostatic  potential  Superimposed  on  each  contour  plot  is  a 

dashed  line  depicting  nQ(y,t)  for  reference  purposes.  Our  initial  electron 
density  profile  nQ(y,0)  is  taken  from  data  supplied  to  us  by  Tsunoda7.  It 
is  derived  from  an  off-perpendicular  VHF  ALTAIR  scan  taken  at  08:04  UT  on  17 
July  1979.  We  found  it  necessary  to  introduce  some  smoothing  of  the  raw  data 
using  a  standard  Shuman  filter  in  order  to  eliminate  spurious  regions  of 
stability  and  instability  in  the  initial  profile.  In  general,  the  bottomside 
gradient  scale  lengths  associated  with  this  profile  are  quite  a  bit  larger 
than  those  we  have  used  in  previous  simulations.4’ 5  This  factor  by  itself 
would  tend  to  give  us  smaller  linear  growth  rates  for  the  collisional 
Rayleigh-Taylor  instability.  However,  this  decrease  in  linear  growth  rate  is 
offset  by  the  fact  that  the  whole  ionosphere  is  somewhat  higher  (this  effect 
tends  to  increase  the  linear  growth  rate)  than  that  used  in  our  previous 
work  (see  reference  4  for  details).  However,  the  higher  altitude  of  the 
ionosphere  will  not  offset  one  other  effect  of  the  large  initial  gradient 
scale  lengths,  viz.,  displacing  a  plasma  fluid  element  vertically  a  given 
distance  will  result  in  a  smaller  relative  depletion  level.  Thus,  all  other 
things  being  equal,  larger  gradient  scale  lengths  in  the  initial  electron 
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density  profile  will  produce  less  depleted  bubbles.  To  summarize,  we  expect 
approximately  equal  growth  times,  but  less  depleted  bubbles,  than  we  have 
seen  in  our  previous  work. 

Our  initial  perturbation  is  given  by  a  pure  sine  wave  of  wavelength 
40  Ax  (our  system  length  in  the  x  direction) : 


Two  simulations  have  been  run:  i)  S,  with  Ax  =  200m;  and  ii)  L,  with  Ax  = 

5  km.  These  two  cases  are  meant  to  span  the  range  of  actual  observed  hori¬ 
zontal  scale  lengths.  Figure  2  shows  isodensity  contours  of  calculation  S 
at  six  times  during  the  simulation.  Figure  3  shows  the  same  contours  at  six 
different  times  for  calculation  L.  The  presence  of  lower  density  plasma  in 
the  bubble  in  calculation  L  is  obvious.  Also  obvious  is  the  fact  that 
calculation  S  seems  to  proceed  in  two  separate  stages,  with  a  small,  low 
depletion  level  bubble  going  through  the  F2  peak  at  about  1200  seconds, 
followed  by  the  main,  somewhat  more  depleted,  bubble  800  seconds  later.  This 
is  in  contrast  to  calculation  L  where  plasma  from  much  lower  altitudes  is 
drawn  up  into  the  bubble  in  one  stage.  Note  that  in  both  cases  we  have  fully 
developed  plumes  (bubbles)  in  2000  seconds. 

Late  time  contours  of  relative  electron  density  n(x,y,t)/nQ(y ,t)  are 
shown  for  calculations  S  and  L  in  Figures  4a  and  4b  respectively.  Solid 
lines  define  depleted  regions,  while  dashed  lines  define  enhancements.  For 
depletions,  the  contour  levels  are  such  that  for  the  first  (outermost)  con¬ 
tour  n/nQ  is  0.5,  and  for  each  succeeding  contour  n/nQ  is  multiplied  by  0.5. 
Thus,  the  third  solid  contour  line  represents  a  value  of  n/nQ  equal  to 


Fig.  2  —  Sequence  of  six  plots  showing  iso-electron  density  contours  of  calculation  S  at  0,  1200, 
I860,  1850,  1950,  and  2050  sec.  Superimposed  on  each  plot  is  a  long  dashed  line  depicting 
nQ(y,t).  Electron  densities  are  given  in  cm'3.  The  observer  is  looking  southward. 
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(e)  (f) 

Fig.  2(Cont’d)  —  Sequence  of  six  plots  showing  iso-electron  density  contours  of  calculation  S  at  0, 
1200,  1650,  1850,  1950,  and  2050  sec.  Superimposed  on  each  plot  is  a  long  dashed  line  depicting 
n0(y.t).  Electron  densities  are  given  in  cm'3.  The  observer  is  looking  southward. 
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Fig.  3  —  Sequence  of  six  plots  showing  iso-electron  density  contours  of  calculation  L  at  0, 1000, 
1500,  1750,  1850,  and  1946  sec.  Superimposed  on  each  plot  is  a  long  dashed  line  depicting 
n0(y,t).  Electron  densities  are  given  in  cm'3.  The  observer  is  looking  southward. 
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Fig.  3(Cont’d)  —  Sequence  of  six  plots  showing  iso-electron  density  contours  of  calculation  L  at  0, 
1000,  1500,  1750,  1850,  and  1946  sec.  Superimposed  on  each  plot  is  a  long  dashed  line  depicting 
n0(y,t).  Electron  densities  are  given  in  cm'3.  The  observer  is  looking  southward. 
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Fig.  4  —  Contours  of  constant  n(x,y,t)/n0(y,t)  for  a)  calculation  S  at  2050  sec  and  b)  calculation 
L  at  1946  sec.  Depletions  (n/nQ  <  1)  are  shown  as  solid  lines  while  enhancements  (n/nQ  >  1)  are 
shown  as  short  dashed  lines.  The  first  (outermost)  depletion  contour  is  for  n/nQ  *  0.5,  while  each 
succeeding  contour  is  for  a  value  of  n/nQ  a  factor  of  0.5  times  the  previous  one.  The  first  enhance¬ 
ment  contour  is  for  n/nc  =  2.0,  while  each  succeeding  contour  is  for  a  value  of  n/nQ  a  factor  of 
2.0  times  the  previous  one.  The  superimposed  long  dashed  line  depicts  n0(y,t). 
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(0.5)3  «  0.125,  or  an  87.5%  depletion.  For  the  first  dashed  contour  line, 
n/nQ  is  2.0,  and  for  each  succeeding  contour  line  n/nQ  is  multiplied  by  2.0. 
Percentage  enhancements  and  depletions  are  obtained  by  subtracting  1.0  from 
n/nQ.  Comparison  of  figures  4a  and  4b  shows  that  the  depletion  level  of  the 
large  horizontal  scale  bubble  is  greater  than  that  of  the  small  horizontal 
scale  bubble.  Furthermore,  neither  of  these  calculations  attain  the  deple¬ 
tion  levels  seen  in  our  previous  work  (see  Figure  7  in  reference  5) .  For 
instance,  in  the  calculation  shown  in  Figure  7b  of  reference  5,  which  is 
identical  to  our  calculation  L  except  for  the  initial  ambient  electron 
density  profile,  a  plume  of  almost  200  km  vertical  extent  can  be  found  with 
depletion  levels  of  99.9%  or  greater.  A  look  at  Figure  4b  shows  that  a  plume 
of  similar  dimensions  can  be  found  only  for  depletion  levels  of  94%  or  greater. 
These  smaller  depletion  levels  for  the  ALTAIR  profile  were  expected  based  on 
our  analysis  of  the  effects  of  larger  gradient  scale  lengths  in  the  initial 
electron  density  profile,  presented  earlier  in  this  section.  An  explanation 
of  why  larger  horizontal  scale  lengths  produce  more  severely  depleted  bubbles 
is  given  in  reference  5.  Briefly,  scale  analysis  is  invoked  to  show  that  the 
vertical  extent  of  the  polarization  electric  field  produced  by  a  perturbation 
in  the  ionosphere  scales  as  the  horizontal  extent  of  that  perturbation. 

Since  it  is  the  electric  field  which  produces  plasma  movement,  the  vertical 
extent  of  the  polarization  electric  field  will  determine  the  depth  in  the 
ionosphere  from  which  a  bubble  may  draw  plasma.  The  lower  in  altitude  that 
a  plasma  fluid  element  originates,  the  lower  its  density  and  hence  the  higher 
the  depletion  level  of  any  plume  into  which  it  is  drawn  (see  reference  5  for 
details).  To  illustrate  this  point,  we  show  in  Figures  5a  and  5b  contours  of 
the  polarization  potential  for  calculation  S  at  1650  seconds  and  for 
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Fig.  5  —  Contours  of  constant  electrostatic  potential  0j  for  a)  calculation  S  at  1650  sec.  and 
b)  calculation  L  at  1750  sec.  Positive  potentials  are  shown  as  solid  lines,  while  negative  potentials 
are  shown  as  short  dashed  lines.  The  contour  levels  are  equally  spaced  from  minimum  to  maximum, 
except  that  the  zero  contour  is  suppressed.  Superimposed  on  each  plot  is  a  long  dashed  line 
depicting  nQ(y,t). 


calculation  L  at  1750  seconds,  respectively.  Contours  of  constant  are 
in  fact  streamlines  for  this  flow  (see  (4)  and  (6)).  Calculation  S  is  seen 
to  consist  of  two  convective  cells,  each  mixing  plasma  over  a  fairly  narrow 
altitude  range,  while  calculation  L  has  formed  a  deep  convective  cell, 
drawing  plasma  into  the  plume  from  very  low  altitudes. 

IV.  SUMMARY 

Both  of  our  global  large  scale  collisional  Rayleigh-Taylor  nonlinear 
simulations,  using  a  background  electron  density  profile  from  ALTAIR  incoherent 
radar  measurements  on  17  July  1979  prior  to  spread  F  onset,  indicate  fully 
developed  ESF  bubbles  in  approximately  one-half  hour.  ALTAIR  radar  and 
ionosonde  measurements  indicate  fully  developed  spread  F  backscatter  about 
one  hour  after  the  profile  was  taken.  The  most  obvious  explanation  is  that 
perhaps  we  have  simply  chosen  our  perturbation  amplitude  too  large;  however, 
since  no  in  situ  measurements  were  made  prior  to  ESF  the  measurements  do  not 
provide  us  with  all  the  pertinent  initial  conditions  prior  to  ESF  onset. 

Other  factors  such  as  the  shorting  effects  of  background  E  region  conductivities 
or  the  neglect  of  inertial  terms  in  the  ion  momentum  equation  could  also 
influence  the  speed  with  which  the  instability  proceeds.  In  addition,  since 
the  chain  of  events  which  leads  from  kilometer  scale  bubbles  to  one  meter 
backscatter  irregularities  is  not  well  understood,  neither  are  the  associated 
time  delays.  Nevertheless,  our  simulations  exhibit  results  which  are  con¬ 
sistent  with  the  onset  time  of  ESF  during  the  17  July  1979  Kwajalein  campaign. 

A  word  of  caution  is  in  order  here  with  regard  to  a  comparison  of  the 
bubbles  we  show  in  our  simulations  and  the  actual  structures  into  which  the 
PLUMEX  I  rocket  was  launched6 * 9  at  12:31  UT  on  17  July  1979.  According  to 
ALTAIR  measurements6,  starting  with  the  onset  of  spread  F  at  about  09:00  UT 


and  lasting  until  approximately  the  time  of  rocket  launch,  the  bottomside  of 
the  F  region  moved  downward  approximately  140  km  at  the  same  time  developing 
much  smaller  gradient  scale  lengths3  than  we  have  used  in  the  present 
simulation.  In  addition,  the  bubbles  detected  by  the  in  situ  plasma  probe  at 
these  later  times  shows  depletions  <  90%.  At  least  three  factors  are  in 
competition  here:  1)  the  downward  movement  of  the  ionospheric  plasma  indicates 
the  presence  of  an  electric  field  which  would  tend  to  reduce  the  growth  rate 
of  the  instability;  2)  lower  altitudes  mean  larger  v  which  again  would 
reduce  the  growth  rate;  and  3)  the  development  of  smaller  bottomside  gradient 
scale  lengths  would  enhance  the  growth  rate.  In  addition,  the  ALTAIR  radar 
operating  in  the  coherent  mode  showed  that  VHF  plumes  (looking  at  lm  field 
aligned  irregularities)  at  this  late  time  during  ESF  were  in  a  decay  phase  of 
development.  Decaying  plumes  which  had  been  generated  early  while  the  iono¬ 
sphere  was  high  would  not  disappear  simply  because  the  plasma  had  been  dis¬ 
placed  downward  later  in  the  evening.  Obviously  the  actual  physical 
situation  at  the  rocket  launch  time  is  more  complex  than  that  addressed  by 
the  simulations  presented  here,  although  many  of  these  complexities  will  be 
treated  accurately  in  forthcoming  versions  of  our  simulation  model. 
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01CT  ATTN  0.  CROMB1E 
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OICY  ATTN  F.  RE  I BERT 
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OICY  ATTN  DAS  I  AC 
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OICY  ATTN  G.  HILLMAN 
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P.  0.  BOX  6770 
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OICY  ATTN  MARSMAL  CROSS 

ILLINOIS,  UNIVERSITY  OF 
DEPARTICNT  OF  ELECTRICAL  ENGINEERING 
URBANA,  IL  61803 
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DEL  HAR,  CA  92014 

01CY  ATTN  S.  R.  GOLDMAN 


JOHNS  HOPKINS  UNIVERSITY 
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JOHNS  HOPKINS  ROAD 
LAUREL,  HD  20810 

01CY  ATTN  OOCUXENT  LIBRARIAN 
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LOCKHEED  MISSILES  £  SPACE  CO  INC 
P.  0.  BOX  50l 
SUNNYVALE,  CA  94088 

01CY  a’-tn  DEPT  60-12 
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P.  0.  BOX  7463 
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OICY  A”N  *.  "EAGHER 

Ll^ACABl’  CORP 
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SAN  DIEGO,  CA  92121 

OICY  ATTN  IRWIN  JACOBS 

M.I.T.  LINCOLN  LABORATORT 
P.  0.  BOX  73 
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OICY  AT-N  DAVID  H.  TOWLE 
OICY  A-TN  P.  WALDRON 
OICY  ATTN  L.  LDUGhLIN 
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5301  BOtSA  AVENUE 
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O’.CT  ATTN  W.  OLSON 
OICY  ATTN  R.  W.  HALPRIN 
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Q1CY  ATTN  P.  FISCHER 
OICT  ATTN  N.  F.  CREVlER 
OICY  ATTN  STEVEN  L.  GUTSCHE 
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PACIFIC-SIERRA  RESEARCH  CORP 
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SANTA  MONICA,  CA  90404 
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few  YORK,  NT  10023 
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01CY  ATTN  £.  A.  STRAKER 
QICT  ATTN  CURTIS  A.  SMITH 
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OICY  ATTN  RAY  L.  LEADA8RAND 
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OICY  ATTN  V.  GONZALES 
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OICY  ATTN  W.  P.  BOQUIST 

TRW  DEFENSE  (  SPACE  SYS  GROUP 
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REOONDO  BEACH,  CA  90278 

OICY  ATTN  R.  <•  PLEBUCH 
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OICY  ATTN  D.  DEE 
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OICY  ATTN  J.  W.  CARPENTER 
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